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With the discovery of the binary black hole coalescence GW150914, the era of gravitational-wave 
astrophysics has started. Gravitational-wave signals emitted by compact binary coalescences will be 
detected in large number by LIGO and Virgo in the coming months and years. Much about compact 
binaries is still uncertain, including some key details about their formation channels. The two 
scenarios which are typically considered, common envelope evolution and dynamical capture, result 
in different distributions for the orientation of the black hole spins. In particular, common envelope 
evolution is expected to be highly efficient in aligning spins with the orbital angular momentum. In 
this paper we simulate catalogs of gravitational-wave signals in which a given fraction of events comes 
from common envelop evolution, and has spins nearly aligned with the orbital angular momentum. 

We show how the fraction of aligned systems can be accurately estimated using Bayesian parameter 
estimation, with 1 cr uncertainties of the order of 10% after 100-200 sources are detected. 

PACS numbers: 04.30.-w,04.80.Nn,04.30.Tv 


INTRODUCTION 

Advanced LIGO [1] has just completed its first science 
run. Results from the first 5 weeks of data have been 
made public, and include the binary black hole (BBH) co¬ 
alescence GW150914 [2]. Advanced Virgo [3] is expected 
to join LIGO in the second science run [4], starting in Fall 
2016. KAGRA [5] and LIGO India [6] are expected to 
join the network before the end of this decade. Ground 
based interferometers will detect gravitational radiation 
from several kinds of sources and start gravitational-wave 
(GW) astrophysics. Gompact binary coalescences (GBG) 
of two neutron stars (BNS), two black holes or a neutron 
star and a black hole (NSBH) have traditionally been 
among the most promising sources, and will be detected 
at a rate of several tens per year, although significant 
uncertainty still exists in the astrophysical rates, both 
predicted and measured [7, 8]. Analysis of many such 
signals promises to shed light on several open problems 
in astrophysics. For example, direct mass measurement 
with GWs could allow for an accurate reconstruction of 
black hole and neutron star mass functions. GWs also 
represent our first chance to perform strong field tests of 
general relativity. An idea of what can be done is given by 
the recent analyses of GW150914 [9-11]. In the last few 
years a pipeline has been created, called TIGER [12-14], 
which can look for unmodeled deviations from general 
relativity in detected GW signals. Elsewhere [15] it has 
been shown how GWs could be used to test proposed 
equations of state for neutron stars. Most of this work, 
and others [16, 17], relies in Bayesian model selection, 
instead of simple parameter estimation. This has the 
advantage that information can be aggregated from all 
detected signals, in a cumulative way, resulting in more 


powerful tests. In this letter we show that GWs can 
be used to check whether spins are preferentially aligned 
with the orbital angular momentum in BBH and NSBH 
systems (we will not consider BNS, since known neutron 
stars in binaries do not have large spins [18]). This is of 
fundamental importance for astrophysics and for under¬ 
standing the formation mechanisms of compact binaries, 
for which there still are many open questions [19-25]. It is 
believed that two main formation channels exist for com¬ 
pact binaries (See [11] for a review). Gommon envelope 
evolution is expect to happen in galactic fields, whereas 
dynamical capture could happen in dense environments 
such as globular clusters. Gritically, it has been sug¬ 
gested that common envelope evolution in binaries will 
align the spins with the orbital angular momentum [26] ^. 
Spins are instead expected to be randomly oriented for 
GBCs formed dynamically. Ultimately, being able to ver¬ 
ify if and how often spins are aligned could significantly 
help to understand the formation patterns of binary sys¬ 
tems, verify which channel happens more frequently, and 
the efficiency of the common envelope evolution phase in 
aligning spins. In this paper we consider a scenario where 
a fraction fa of signals have spins nearly aligned, while 
the rest are non-aligned. This accounts for two possible 
formation patterns for GBCs, one of which is efficient in 
aligning spins with the orbital angular momentum. We 
find that the posterior distribution for the mixture pa¬ 
rameter fa can be accurately estimated with a couple 


^ Others suggest that kicks introduced in the system when the 
progenitor stars undergo core collapse supernovae could result in 
spins being significantly misaligned [27] 



2 


hundred sources, with precision of the order of 10%. 


Modulo a normalization constant, the log of Eq. 2 then 
reads: 


METHOD 

Let us assume two main formation channels for CBC 
exist, which result in a fraction fa of systems having spins 
nearly aligned, and a fraction {1-fa) having misaligned 
spins. We will assume N GW detections are made, de¬ 
noted by their data streams di with i = 1 ... N, and show 
how they can be used to estimate fa- 

We introduce two mutually exclusive models: Ha cor¬ 
responding to spins nearly aligned with the orbital an¬ 
gular momentum, and Ha corresponding to non-aligned 
spins (we will define these models more precisely later). 
Given an event and the corresponding data stream dk, 
we can calculate the evidence of the data for the mod¬ 
els above, ZJf = p{dk\Hm), with m = a,d. The evidence 
must be calculated by integrating a non-trivial likelihood 
function over a multidimensional parameter space [28]. 
Galling 9 the unknown parameters on which a GBC de¬ 
pends, we have: 


Z]f= / p{dk\eHm)Pm{9\H^)de . (1) 

We solve this integral by using the Nested Sampling 
and the BAMBI flavors of lalinference [28]. We stress 
that the hypervolume we integrate over, 0^, depends on 
the model being considered (see next section). 

We can now show how the aligned fraction of events 
can be calculated. We start by applying Bayes’ theorem 
and the product rule to the posterior distribution of fa- 

N 

p{fa\d) (X p{d\fa)p{fa) = p{fa) Wp{dk\fa) (2) 

where p{fa) is the prior on the mixture parameter, that 
we take as flat on [0,1] since we do not have any previous 
astrophysical information. With d we have denoted the 
set of N detected events, d = {di, ^2, ■ ■ • , dAr}- The fac¬ 
tors inside the product can be expanded by noticing that 
Ha and Ha are mutually exclusive for each detection, i.e. 
that p{Ha) -\-p{Ha) = 1. Thus: 


p{dk\fa)= p{dk\Hjfa)p{Hj\fa) (3) 

j=a,a 

We notice that Hjfa can be written simply as Hj, 
since knowing if a signal was aligned or not makes know¬ 
ing fa irrelevant. Next, we need to calculate p{Ha\fa) 
and p{Ha\fa)- These are trivially p{Ha\fa) = fa and 
p{Ha\fa) = {1 — fa)'- if & fraction fa of events is aligned, 
the probability that the aligned model applies to any 
event, before looking at the data, is fa- 


logp{fa\d) 


log(p(/a)) (4) 

N , 

^(logZ,“ + log 

fc=l 



IMPLEMENTATION 

In this section we describe the parameters of the 
sources we simulate. We consider GWs emitted by BBH 
and NSBH, and for each type of source we generate 
two catalogs of GW signals, one corresponding to nearly 
aligned spins and the other to non-aligned spins. 

For the BBH, we use the so called IMRphenomPv2 
waveform approximant [29, 30] (this is one of the two 
families used for the analysis of GW150914 [9]). For 
the lighter NSBH we use the SpinTaylorTd approxi¬ 
mant [31, 32]. Unlike IMRphenomPv2, SpinTaylorTd is 
an inspiral-only approximant and thus cannot model the 
merger and ringdown phase of GBCs. However, since 
the frequency at which the merger happens is roughly 
inversely proportional to the total mass of the system, 
merger and ringdown can be neglected as long as the to¬ 
tal mass is below ~ 20Mq [33], for reasonable signal-to- 
noise ratios (SNRs). All the NSBH we simulate have total 
mass below 13Af©. In both cases, we work at the high¬ 
est known post-Newtonian phase order, while neglecting 
higher-order amplitude corrections. We also neglect tidal 
contributions in neutron stars [15, 34-36]. Both these 
limitations are due to computational considerations and 
will not impact our main result. 

For the BBH, we choose to consider heavy black holes 
of a few tens of solar masses, which we know will be 
detected in large number by ground based detectors 
in the coming months [7]. We thus generate compo¬ 
nent masses uniformly from the range [30, 50]Mq (in 

the source frame). The dimensionless spins, 

are uniformly generated in the range [0,0.98], compati¬ 
ble with the range of validity of the waveform approx¬ 
imant [29, 30]. For the NSBH signals, BH masses are 
in the range [6,11]M0, and NS masses in the range 
[1.2, 2]M©. Black hole spins are uniform in the range 
[0,1] while for the neutron stars we restrict possible spins 
to [0, 0.3] (the largest measured spin of a NS in a bi¬ 
nary system is 0.02 [37]). For the NSBH sources, our 
choice for the mass range of BH is mainly driven by the 
use of inspiral-only waveforms (IMRphenomPv2 wave¬ 
forms are not considered reliable for mass ratios above 
4-5 [29, 30]). 

The distances are uniform in comoving volume, with 
a lower network SNR (that is, the root-sum of squares 
of the SNR in each instrument) cut at 8-\/3 ^ 13.9 for 
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NSBH and 7-\/3 ~ 12.1 for BBH. These correspond to 
distances up to ~ 1.2 Gpc for NSBH and ~ 12 Gpc for 
BBH. For both BBH and NSBH, the sky position and 
orientation-polarization of the systems are uniform on 
the unit sphere. 

To verify that the test we propose is self-consistent 
and does not rely on the exact definition of “aligned”, 
we define it in a different way for NSBH and BBH. 

For NSBH, the nearly aligned (henceforth just aligned) 
catalog is made of signals with tilt angles (i.e. the an¬ 
gles between spin vectors and orbital angular momen¬ 
tum) isotropic in the interval [0,10]°, i.e. close to the 
positive direction of the orbital angular momentum. For 
the BBH, the tilts in the aligned catalog are in the range 
[0,10]° U [170 — 180]°, i.e. the spin vectors can be along 
both the positive and negative direction of the orbital 
angular momentum. For both BBH and NSBH, the non- 
aligned model is the logical negation of the correspond¬ 
ing aligned model. For example, for NSBH tilts were 
isotropic in the range [10,180]°. The priors on the tilt 
angles for the Ha and Ha models, eq. 1, are isotropic 
with cuts that match these intervals. 

Each event is added into simulated Gaussian noise cor¬ 
responding to the design sensitivity of the two Advanced 
LIGO detectors and the Advanced Virgo detector [4]. 

We analyze all events in the catalogs twice, once with a 
prior that matches the Ha model, and once with a prior 
that matches Ha- These runs provide the evidences of 
eq. 1 that we can use in eq. 4 to calculate posterior of the 
mixture fraction. 

RESULTS 

To show how the method performs for some representa¬ 
tive values of /a, we generate (for both BBH and NSBH) 
5 catalogs with increasing fraction of aligned events. 
From 0 (all events are not-aligned) to 1 (all events are 
aligned), with steps of 0.25. These catalogs are trivially 
created from our initial set of NSBH and BBH by ran¬ 
domly drawing aligned and non-aligned signals with the 
desired ratio until 100 sources for NSBH or 200 for BBH 
are obtained. The evidences of these events are then used 
in eq. 4 to obtain the posterior distribution of fa- The 
main results are shown in Fig. 1, where for pedagogical 
purposes we keep separated BBH and NSBH sources. We 
see how the posterior distributions for fa peak at or very 
close to the corresponding true values, given in the leg¬ 
end, with 1 a uncertainties of the order of 10%. The small 
offsets of some of the curves can be explained with the 
limited number of events we consider (the offset would 
be zero in the limit of an infinite number of sources). By 
regenerating the catalogs a few times we saw that the 
peaks can shift by a few percents on either side of the 
true values. We have verified that halving the number of 
sources in the catalogs (50 NSBH and 100 BBH) broadens 


the posterior distributions, while leaving them centered 
around or close to the true values. One might be sur¬ 
prised that the NSBH distributions are narrower in spite 
of the fact that fewer (100) sources are used for NSBH 
than for BBH (200). This can be explained by remind¬ 
ing that for BBH we used a slightly lower SNR threshold 
(12 vs 13.9), thus increasing the number of weak signals 
that do not contribute much to the measurement, while 
broadening the posterior distributions. Furthermore, the 
characteristic effects of misaligned spins (e.g. amplitude 
precession) are more visible in NSBH, which make them 
ideal sources for this test. 



Figure 1: Posterior distribution for the mixture 
parameter fa after 100 NSBH (dashed) and 200 BBH 
(solid) detections. Several underlying values of fa 
(given in the legend) are considered, /a = 0 corresponds 
to a catalog where none of the sources had aligned 
spins, while /a = 1 refers to a catalog where all events 
had aligned spins. 

We stress that we do not assume that the priors in 
eq. 1 perfectly match the corresponding distributions in 
the simulated events, which will likely happen in the 
first years of gravitational-wave astrophysics. Two im¬ 
portant examples are the prior distributions for distance 
and masses. 

While geometrical arguments led us to use a prior for 
the luminosity distance uniform in comoving volume, in 
reality, since far away sources would not be detectable, 
the distance distribution of detected events will first in¬ 
crease with distance, reach a maximum, and then de¬ 
crease. The distances corresponding to the maximum of 
the distribution and the length of the tail depend on the 
true astrophysical distribution of masses (heavy CBG will 
be visible farther away), which we don’t know (but will 
hopefully measure in the coming years). Similarly for 
the mass prior: we used priors in the component masses 
which were a factor of few larger than the range used to 
simulate the sources. It will be interesting to verify how 
the test performs if the true distribution of tilt angles for 
the aligned model is different than what is used to split 
the two models, or if the true distributions are not mu- 











4 


tually exclusive^. Given the amount of simulations that 
would be necessary to fully explore those scenarios, we 
leave it for future work. 

It is worth remarking that while we consider two possi¬ 
ble formation channels, this framework can be extended 
to take into account more models, provided they are mu¬ 
tually exclusive. Similarly, if one believes only one for¬ 
mation channel is possible in the universe, and thus all 
events will either have aligned or not-aligned spins (that 
would correspond to fa = 0,1), model selection can be 
used to quantify how many detections are required before 
the model can be proven right. 

Although we believe considering a mixture of two mod¬ 
els is more consistent with today’s understanding of bi¬ 
naries’ formation, we give an example of a single-channel 
test. For this, we simulate a situation in which all sources 
have non-aligned spins and we calculate the cumulative 
odds ratio: 

Qa ^ Pj'Hald) ^ p{d\na)p{na) ^ pCHg) 
piV-ald) pid\'Ha)p{'Ha) “ Pi^da) 

where S? is the cumulative Bayes factor for aligned vs 
non-aligned models. Since the data corresponding to the 
N detections is statistically independent, the cumulative 
Bayes factor can be written as a product over the single 
events: 


^ TT PidklTda) 


N 

n 


Z^k 


( 5 ) 


The logarithm of the odds ratio is shown in Fig. 2 as 
a function of the number of events, for random sub¬ 
catalogs of 10 NSBH and 50 BBH. We have assumed 
p{Tda) = piTda)- We see that for both type of sources the 
correct non-aligned model is favored in a significant way 
(log odds below the solid horizontal line favor the non- 
aligned model at a > 2.7a level). NSBH curves go down 
faster since for NSBH the effect of spin misalignment are 
stronger in the waveform, and thus harder to match with 
an aligned model. 


CONCLUSIONS 

Two formation channels are commonly considered for 
CBC: common envelope evolution, which should result in 
spins to be preferentially along (or very close to) the di¬ 
rection of the orbital angular momentum, and dynamical 


^ In our simulations we assumed isotropic tilt distributions p{t) oc 
sinr, with cuts at 10°. This makes our non-aligned distribu¬ 
tion basically equal to a fully isotropic distributions, since the 
probability that both tilts are small (or close to tt for BBH) is 
negligible. 
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Figure 2: (color online) Cumulative odds ratio for 
NSBH (red) and BBH (blue), with non-aligned 
injections. Each line is a sub-catalog. Cumulative odds 
values below the solid horizontal thick line favor the 
(correct) non-aligned model with a significance larger 
than ~ 2.7 a. We notice that we cut the y axis at -250 
to improve clarity. Some NSBH catalogs go down to 
cumulative odds of -1000. 


capture, which should results in randomly oriented spins. 
In fact, there is not complete agreement on whether com¬ 
mon envelope evolution is efficient enough in aligning 
spins, or if instead eventual kicks from the core collapse 
supernova of the progenitor stars will be the dominant 
factor. It would thus be of importance to calculate which 
fraction of the compact binaries have spins nearly aligned 
with the orbital angular momentum, which could be used 
to expand our understanding of formation channels. In 
this paper, we have shown how gravitational waves emit¬ 
ted by compact binaries containing a black hole, and de¬ 
tected by Advanced LIGO and Virgo, can be used to 
verify if spins are preferentially aligned with the orbital 
angular momentum. We considered neutron star - black 
hole and binary black hole systems, and created cata¬ 
logs of sources with increasingly large fraction of aligned 
sources (from 0 to 100%). Black holes in NSBH were 
of low mass (up to IIMq), while for BBH we simulated 
heavy objects, comparable to GW150914 (masses in the 
range [30 — SOJM©), which will be detected in large num¬ 
ber in the coming months and years. 

We showed how a couple hundred signals are enough 
to pinpoint the underlying value of the aligned fraction 
with ~ 10% uncertainty, which suggests GWs represent 
a viable way of gaining insight into the orientation of 
spins in compact binaries, and ultimately on their evolu¬ 
tion. We have verified the robustness of the test against 
some common prior mismatch (distance, masses). Future 
work includes introducing a mismatch between the defi¬ 
nition of aligned in the test and the true distribution of 
aligned sources. We also stress that if more information 
is available which could help distinguish between the two 
channels (e.g. the resulting mass ratio distribution), it 
could be folded in an extended version of this test. 
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